Measurements of Intensity Distributions in the Approach to Localization 



M. Stoytchev and A. Z. Genack 
Dept. of Physics, Queens College of CUNY, 65-30 Kissena Blvd., Flushing, NY 11367 

(April 10, 1998) 

We report measurements of intensity distributions of transmitted microwave radiation 
in quasi- ID samples with lengths L as large as the localization length £. In contrast to neg- 
ative exponential statistics found in the diffusive limit, the distribution falls as a stretched 
exponential of power 1/2 and becomes nearly log-normal as L approaches £. We confirm 
the relation between the moments and full distributions of the intensity and total transmis- 
sion obtained by Kogan and Kaveh using random matrix theory. Good agreement is found 
when this relationship is used to compare measurements of intensity distribution in strongly 
absorbing samples with lengths as large as £ with calculations by Brouwer of the transmis- 
sion distribution. The variances of the measured distributions increase superlinearly near 
the localization threshold confirming the coexistence of localization and absorption and 
the suitability of these variances as measures of the closeness to the localization threshold. 
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As the localization threshold is approached, fluctuations in key transmittance quantities become as large 
as the ensemble average values of these quantities. Thus a comprehensive description of mesoscopic trans- 
port should provide the full distribution of transmittance quantities and the relationships between them. |l| 
In order of increasing spatial averaging, key transmittance quantities are the intensity, corresponding to the 
transmission coefficient T ao for incident mode a and outgoing mode 6, the total transmission, T a — ^2 b T a {>, 
and the transmittance, T = ^2 ab T ao . In the absence of absorption, the transmittance equals the dimen- 
sionless conductance g = Nl/L, where ./V is the number of transverse modes and I is the transport mean 
free path, g] Localization in quasi-lD samples is achieved when <7 = latL = £ = Nl. In this Letter, we 
focus on the intensity distribution, which is the key distribution in statistical optics. Q| We demonstrate 
its relationship to the distribution of total transmission, find the scaling of the variance of the intensity 
and total transmission up to L — £, and determine the extent to which absorption influences localization. 

In the diffusive limit, the degree of long-range intensity correlation is small and the intensity distribution 
is well approximated by the Rayleigh distribution. [||-|5| For polarized detection, this corresponds to negative 
exponential statistics, P{s ab ) = exp(— s ab ), where s ab = T ab / < T ab > is the intensity normalized to 
its ensemble average value. In previous work, deviations from negative exponential behavior have been 
observed and ascribed to long-range intensity correlation ^,^,0. In these studies, fluctuations as large as 
s a b ~ 10 ~ g were observed. In the present work, fluctuations as large as fifty times g are observed in 
samples with lengths L ~ £. 

The intensity distribution is studied in a quasi-lD geometry, which is equivalent to the electronic case 
of a thin wire. Thouless argued that, the level width 5v in a wire at T = should become smaller than 
the level spacing Av since 5v is proportional to the inverse of the travel time and so falls as 1/L 2 in the 
diffusive limit, whereas Av is the inverse of the density of states and so falls as 1/L. As a result, the modes 
in adjacent sections of the wire should not overlap, and electrons should become localized. || 

The question arises as to whether radiation can be localized in the presence of absorption. In this 
case, the level width falls as 1/L, just as the level spacing does. In previous measurements in absorbing 
samples, the variance of the normalized total transmission, < s ab >= T a / < T a >, var(s a ), was found 
to scale sublinearly with L. || If the attenuation length due to absorption, L a , serves as a cutoff length 
for localization p~0| . |TT[ | then var(s a ) would approach an asymptotic limit as L increases. On the other 
hand, if localization can be achieved in absorbing samples, then var(s a ), which is essentially the degree of 
correlation in the intensity of different outgoing modes, should increase superlinearly as L approaches £. 
This might occur, since the wave remains temporally coherent in the presence of absorption. Jl2|,[l3| Weaver 
has shown in a 2-D simulation that the introduction of absorption does not disrupt the spatial localization 
of acoustic waves in closed systems, though the overall energy decreases exponentially with time. In 
recent calculations, Brouwer found that for diffusive waves the prefactor multiplying L/£ in the expression 
for var(s a ) drops from | to ^ as the ratio L/L a increases. [ fl4| The behavior of this quantity, however, was 
only considered for lengths considerably less than the localization length. 
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Here we report measurements of intensity transmitted through random waveguides with L < £, but 
> L a . We expect that modes in this sample are completely mixed and the degree of intensity correlation 
between different modes is constant. Wave propagation in this sample should therefore be described by 
random matrix theory (RMT). jL5| Recently, Kogan and Kaveh used RMT to obtain a relationship between 
the moments of intensity and total transmission in nonabsorbing quasi-lD samples. ]l6[ ] They find, 

< S ' Q l b >=n!< S >, (1) 

This leads to a relationship between the distributions of intensity and total transmission, |l6| , |T7[ 

P( Sab ) = / ^ P ( Sa )exp(- S -^). (2) 

Since the distribution of total transmission can be calculated from the distribution of the eigenvalues of 
the transmission matrix, these relations provide a basis for calculating the intensity moments 

and distributions from RMT. 

The intensity distribution is obtained from measurements of the field transmitted through an ensemble 
of random configurations of 1.27 cm-diam. polystyrene spheres inside a copper tube. Tubes with diameters 
d = 5.0 and 7.5 cm and lengths up to 520 cm are used. The samples have filling fractions of 0.52 and 
0.55 for d = 5.0 and 7.5 cm, respectively. The sample tube is rotated between successive measurements to 
produce new configurations of scatterers. At least two thousand sample configurations were used for each 
distribution. Field spectra were taken from 16.8 to 17.8 GHz in steps of 1 MHz using a Hewlett-Packard 
8722C network analyzer. The radiation is coupled into and out of the sample by 0.4 cm wire antennas 
placed 0.5 cm from the ends of the sample. In order to ensure that the distributions were not distorted by 
noise, it was necessary to use an amplifier with an output power of 40 W for samples with lengths greater 
than 200 cm so that the average intensity was at least three hundred times the noise. In the frequency 
range of the measurements, I ~ 5 cm |2Cj. A fit of measurements of the field autocorrelation function 
with frequency shift to theory Q gives L a — 34 i 2 cm and D = (3.03 ± 0.21) x 10 10 cm 2 /s, where D 
is the diffusion coefficient. The localization length for the samples with d = 5 cm is found below to be 
£ = 551 ± 18 cm. 

In Fig. 1, we present the intensity distributions for two samples with L/£ ~ 0.1 and 0.4. Calculations 
for diffusive waves in the absence of absorption have predicted that for s a b 3> Sj the intensity distribution 
falls as exp(—2y/gs a b). (l8) For the samples measured, we find P(s a b) ~ exp(—2y/js a b) in the tail of the 
intensity distribution. The values of 7 obtained from a fit to the tail of the measured distributions is within 
20% of the parameter g' = 3va ^ s ) , || which equals g in the absence of absorption . The fit of a stretched 
exponential to the distribution for L/£ = 0.4 in the range of s a b from 10 to 18 is shown in Fig. 2 and gives 
7 = 2.9, which is close to value of g' of 3.06 for this sample. 

The measured intensity distributions are compared to the transform of the measured transmission distri- 
butions H for the corresponding samples using Eq. (0). The transforms shown as solid lines in Fig. 1, are 
in good agreement with the measured intensity distributions. A comparison of the moments of intensity 
and transmission is shown in Fig. 3. We find an increasing deviation of the ratios <s" b > /n! <s"> from 
the value of unity expected from Eq. ([!]) as n increases. We find, however, that agreement with Eq. (]l|) is 
dramatically improved when the moments are calculated using the asymptotic expressions in the diffusive 
limit for the intensity and transmission distributions beginning from the point at which the measured 
distribution have their first zero. The asymptotic expressions for the intensity distribution exp(—2^/g's a ii) 
is substituted for the measured distribution for values of s a b between 20 and 150, whereas the asymptotic 
exponential expression exp{—g's a ) P Jl^ , [l^| is substituted for the measured transmission distribution for 
values of s a between 5 and 25. The improved agreement indicates that the extent to which the measured 
ratio of moments is in accord with Eq. (Q) is largely limited by the range of intensity and transmission 
values measured, which depends on the number of configurations on which measurements were made. 

Applying Eq. (Q) to the second moments gives 

var(sab) = 2var(s a ) + 1. (3) 



The same expression can be obtained from a perturbation calculation up to order 1/g [ jii||22] ]. Our mea- 
surements confirm the prediction of RMT that Eq. ([|) is independent of the value of g and is correct up 
to order 1/N. 
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Intensity statistics at the localization threshold are studied in measurements at L = 520 cm and d = 5 
cm and shown in Fig. 4a. Values of s a b as large as 50 are observed. The distribution is seen in Fig. 4b to 
be nearly log-normal, in agreement with predictions for localized radiation JUIO! ■ 

Using Eq. (^) , we compare these measurements of the intensity distribution to random matrix calculations 
of the transmission distribution in the presence of absorption G3]. In order to compare the intensity 
distribution to theory, however, £ must be determined. Far from the localization threshold, in the absence 
of absorption and internal reflection, var(s a ) and £ are related by Jl4] , ^6| , p^p^ ] 

var(s a ) = |-. (4) 

To find £ in samples in which corrections due to absorption, localization and internal reflection cannot be 
ignored, we first obtain the intensity distribution for the equivalent samples without absorption using the 
measured spectra in our absorbing samples. The procedure used is based on the work by Yossefin ]l3| who 
proposed that for ujr a 3> 1 and when the dispersion is neglected the field in absorbing media differs from 
that in the absence of absorption only by a factor of exp(t/2r a ), where r a = L a 2 /D s is the exponential 
attenuation time due to absorption. We first obtain the time response E(t) to a narrow gaussian pulse in 
time by Fourier transforming the product of the measured spectrum and a broad gaussian in the frequency 
domain. The time dependent field is then multiplied by exp(t/2r a ) and the so modified time spectra are 
transformed back into the frequency domain. The intensity distribution is then calculated using these field 
spectra. This procedure could be applied, however, only for L < 150. For longer samples the signal in time 
is buried in the noise level for large t and thus a significant part of the field in the absence of absorption 
cannot be recovered. The intensity distributions corrected for absorption are in good agreement with 
transforms of the diffusive result for the distributions of total transmission calculated in Ref. |l6|,[l8| and 
give values for the parameter g' |1 equal to the value of g, as expected in the absence of absorption. We 
also find that the average transmission obtained from the spectra corrected for absorption is consistent with 
the expected scaling as (L + 2zb)~ 1 , where Zb is the diffusion extrapolation length due to internal reflection 
p3| . These results confirm the ability of this approach to statistically eliminate absorption in the diffusive 
limit. The influence of internal reflection upon var(s a ) in the absence of absorption can be accounted for 
by substituting L = L + 2zb for L in Eq. (||). We next account for the leading order correction to var(s a ) 
due to nonlocal correlation. In the absence of absorption, the variance is increased by an additional factor 
of (1 + ||) ||J2| to yield, 

, , 2L 4£ 2 
Var{Sa) = 3? + W (5) 

A fit of Eq. (|) to the data corrected for absorption using £ and Zb as fitting parameters gives £ = 551 ± 18 
cm and Zb — 5.25 ± 0.31 cm. The value of Zb obtained is consistent with the values of this parameter for 
the same samples in the frequency range between 18 and 19 GHz. [f2of 

The dependence of var(s a ) upon L with and without absorption is shown in Fig. 5. The solid curve 
represents the result of the calculations in the diffusive regime (L/£ <C 1) by Brouwer [[lij which account 
for absorption and the thin dashed curve shows the fit of Eq. (^) to the data corrected for absorption. The 
values of var(s a ) are calculated from var(s a b) using Eq. (||). For lengths up to 200 cm, the result from 
the measurements, show sublinear behavior which is consistent with the results from total transmission 
measurements || and the calculations in Ref. ^M- The deviation from the solid line increases for larger 
lengths and may reflect localization corrections that were not included in the theory. For strongly absorbing 
samples (L 3> L a ), Brouwer finds a log-normal distribution for the total transmission with < lnT a > = 
—L/L a — 3L/4^ — lnN and var(lnT a ) — L/2£. (l4j Using thiese relations, we calculate the values of var(s a ) 
in the regime of strong absorbtion. The result is shown by a thick dashed line. We note that this curve 
has a physical meaning for L 3> L a only. In this regime, it shows significant corrections due to localization 
in qualitative agreement with the results from the measurements. Thus, we can associate the increase of 
£var(s Q ) / 'L for these samples with the transition to localization. 

We now compute P(s a b) for the sample with L = 520 cm using the transmission distribution calculated 



in Ref. |14| and the values of £ and L a found here. The calculated intensity distribution is presented as 
the solid line in Fig. 4a and is in good agreement with the measurements. 

In conclusion, we find that the intensity distribution for s a b S> g is described by a stretched exponential 
to power 1/2 and that the distribution for L ~ £ is close to log-normal. We confirm experimentally 
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the relationships obtained by Kogan and Kaveh between the moments and full distributions of intensity 
and total transmission. These relations unify the statistical description of local and spatially averaged 
transmittance quantities. Our measurements demonstrate that the statistics of wave transport is only 
marginally affected by absorption and that absorption does not substantially inhibit localization. The 
ability to reach the localization threshold using a quasi-one-dimensional sample is an extension to classical 
waves of the suggestion by Thouless that electrons will always be localized in sufficiently long wires at low 
temperatures. These results show that the variances of the intensity or transmission are reliable measures 
of the impact of localization upon transport in random media. 
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FIGURES: 



Fig. 1. Probability distribution of intensity for samples with £/£ ~ 0.1 and 0.4; the samples dimensions 
are (a) d = 7.5 cm, L — 100 cm, and (b) d = 5.0 cm, L = 200 cm, respectively. The solid lines represent 
distributions obtained from measured transmission distributions [9] using Eq. (2). 

Fig. 2. Fit of a stretched exponential of power 1/2 to the tail of the intensity distribution (L/£ ~ 0.4). 

Fig. 3. Comparison between the moments of intensity and total transmission ~ 0.4): • moments 
obtained from the measurements, o moments calculated from the extended distributions. 

Fig. 4. Intensity distribution for the sample with L = 520 cm and d — 5.0 cm. The solid line in part (a) 
shows the distribution obtained from a transform of the total transmission distribution for this sample cal- 
culated using the expressions from Ref. [14]. The dashed line in part (b) represents a normal distribution. 

Fig. 5. Dependence of var(s a ) upon L. The different symbols represent: • results obtained from the 
measurements, o results from the data corrected for absorption. Note that the curves obtained using the 
theoretical expressions in Ref. [14] represent the ratio £var(s a )/L, not £var(s a )/L. 
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